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Figure 1. Ground state effective bare binding en- 
ergy. (Locai sink, smeared source). 
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Figure 2. IS, 2S, and 3S wave functions, normal- 
ized to unity at tlie origin. 



Figure 1 displays an effective-f*^ plot for the 
local-smeared correlation matrix. There is a long 
plateau beginning at very early times indicating 
the effectiveness of the procedure in isolating the 
ground state. Figure 2 displays Bethe-Salpeter 
amplitudes for the lowest three S'-wave states, 
normalized to unity at r = 0, and measured us- 
ing dynamically-improved smeared sources. To 
obtain values for the decay constant, correlated 
fits in t and extrapolations in k to Kg = 0.157 
were made with covariance matrices estimated 
from single-elimination jackknives (nested when 
required). Quoted errors are for the bare lattice 
quantities only, are statistical only, and were ob- 
tained by bootstrapping, with 250 samples, the 
entire fitting procedure. We find provisional val- 
ues at 13 = 6.0 of ZlV^ = 0.279te, 0.244^^, 
and 0.224t^ at k = 0.152, 0.154, and 0.155 re- 
spectively. With the renormalization factor [5] 
^stat = Z^2''^V2^^/4:Kc/k - 3 with Z^'^' = 0.70, 
ruB = 5.28Gev and a nominal value of a~^ = 
2.0GeV [6], the extrapolation to k = Kg = 0.157 



fstat 
JB 



2.0GeV 



3/2 



224l^(statistical) (12) 



This is in agreement with the results of the Fer- 
milab group who also use nearly perfect sources, 
but is lower than the results of other groups [7]. 



In conclusion, the variational basis of MOST 
generates optimal interpolating fields objectively 
from the data - no ansatz is needed. These 
are used to determine wave functions, decay con- 
stants and bare binding energies £° . Effective-5'^ 
plots plateau early, are fiat and have small statis- 
tical errors. Results have been obtained for the 
lowest S'-wave states; the method generalizes for 
orbital excitations. The interpolating fields will 
be optimal choices for use in matrix element cal- 
culations of further properties of 5 mesons. 
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entire correlation matrix. This is a unique luxury 
afforded in the static approximation and provides 
the maximal (valence quark) basis to exploit vari- 
ational methods. 

It is illuminating to use matrix notation, and 
to write the correlation matrix as 



C(t) 



Ae 



-Mt 



A^ 



(4) 



with Ain = (0 I $j(0,0) I n). In particular, An 
is the coupling of the local current to the ground 
state n = I. The pseudoscalar decay constant is 
determined from 



Zr 



An = fT'V^^^iZT'r^a'"'' 



(5) 



Zf\ defined by A™"' = Zf^'A^'^'a-^l'^ is the 
renormalization factor [5]. (The^on Z^ and Z^'^' 
indicates that all K-dependent normalizations are 
absorbed into the renormalization factor.) 

However, the ^'s are not orthogonal, and can- 
not be determined directly upon diagonalizing the 
matrix at time slice to >Q 



C(to) = V(to)A(to)VHto) 



(6) 



where A(to) is the diagonal matrix of eigenvalues. 
These are approximately equal to the eigenvalues 
of the transfer matrix; corrections are damped 
exponentially only after large to by which time 
the signal-to-noise ratio has been degraded. 

Kronfeld [3] and Liischer and Wolff [4] have 
outlined procedures by which eigenstates of the 
transfer matrix can be estimated at early times; 
the formalism, adapted to our variational basis 
of all relative spatial separations of the quarks, is 
very effective in isolating eigenstates. Following 
Kronfeld, there exists an orthogonal matrix R(to) 
which relates the eigenvalues A„n(to) and eigen- 
vectors Vin(to) of the correlation matrix to those 
of the transfer matrix: 

(^e-^*°/2) = {vito)Aito)+'^'') Rito) (7) 

R(to) can be determined from the correlation ma- 
trix at later times t > to > 0- Diagonalizing 

D(t,to) = A(toy^V(toyC(t)V(to)A(toy^ 
= i?(to)e-^(*-*°)i?(to)t (8) 

determines the amplitudes 
A = V{to)A{to)+^I^R{to)e+^'''l^ (9) 



from which the decay constant(s) can be calcu- 
lated directly using eqn. 5, and determines Ainn 
(which for static-light mesons are the bare bind- 
ing energies 5^). 

As a important bonus, one obtains 

Z = V{to)A{to)-^I^R{to)e-^'''l^ (10) 

which projects out eigenstates of the transfer ma- 
trix: Z'^ A = 1. Thus $' = Zt$ are optimally im- 
proved interpolating fields, obtained dynamically. 



Z^C{t)Z 



-Mt 



(11) 



We have suppressed, for clarity, factors of n,-, 
the number of spatial points equivalent to a given 
point ri under the cubic symmetry group; these 
remain upon condensing the full L^ x L^ correla- 
tion matrix down to size N x N . 

All two-quark interpolating fields do not pro- 
vide a complete basis (due to multiquark and hy- 
brid intermediate states), so the exact expression 
in eqn. 8 is corrected by exponentially damped 
terms, which however are expected to be smaller 
than those correcting the extraction of eigenval- 
ues of the transfer matrix directly from the eigen- 
values of the correlation matrix [4]. Thus for 
moderate values of t and to one expects eqns. 9 
and 10 to be good approximations. 

The smaller eigenvalues of the correlation ma- 
trix are noisy and contribute little to the largest 
eigenvalues of the transfer matrix. We truncate 
the correlation matrix to its highest n\ eigenval- 
ues before the construction in eqn. 8 [3]. This 
introduces further exponentially damped correc- 
tions which are expected to be small and which 
are dealt with as above. 

For this report, rather than use eqn. 9 di- 
rectly to determine f§^*', we instead use the 
dynamically-improved interpolating fields $' = 
Zt$ (obtained from diagonalizing _D(3,2) with 
nx = b), in the conventional way to determine Z^ 
from fits to local-smeared and smeared-smeared 
correlations. 

Quark propagators were computed for k = 
0.152, 0.154, and 0.155 with the Wilson fermionic 
action in a background of 32 quenched (Coulomb 
gauge-fixed) configurations (selected every 1000 
pseudo-heat-bath sweeps after 5000 thermaliza- 
tion sweeps) with /3 = 6.0 on a 20^ x 30 lattice. 
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We present a new method for the study of heavy-light mesons in the static approximation of lattice QCD 
which is optimally effective in isolating ground and excited states. With "MOST" (Maximal Operator Smearing 
Technique), the heavy quark is smeared at all possible positions relative to the light quark, subject to the constraint 
of cubic symmetry. With correlation functions constructed using this set as a variational basis, eigenstates of the 
transfer matrix are projected out at very small time separations, where statistical errors are small. We illustrate 
the utility of the method with preliminary results for the meson decay constant fg^^, binding energies and wave 
functions of the lowest states. The method produces dynamically-improved interpolating fields which can be used 
for matrix element calculations. 
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Lattice QCD can make important contribu- 
tions to phenomenology by the reliable calcula- 
tion of decay constants, mixing amplitudes and 
form factors of 5-mesons. It is natural to study 
the effective theory where the heavy h quark is 
approximated by a static color source [1] giving a 
simpler system with only one dynamical (valence) 
quark. 

To satisfactorily extract hadronic ground state 
properties from lattice correlation functions of 
operators, it is imperative to choose operators 
which couple much more strongly to the ground 
state than to excited states. With the first at- 
tempts at a lattice calculation of fg^*' in the static 
approximation, it was realized that local opera- 
tors were very poor choices. A second generation 
of attempts used operators which were spatially 
smeared; this had been done successfully for light- 
light mesons. Smearing over cubes, over walls, 
or with Gaussian or exponential functions im- 
proved matters by improving statistics, but were 
still unsatisfactory [1]. The exponential degrada- 
tion of the signal-to-noise ratio with increasing 
Euclidean time, being much more pronounced for 
static-light than for light-light mesons, requires 
that better, more physical, smearing functions be 
used. 
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The use of a spinless relativistic quark model 
Hamiltonian to construct a basis set of trial 
smearing functions for a variational calculation 
was an important improvement [2]. We propose 
an alternative: that the variational calculation be 
done in a manner free of any such ansatz, allow- 
ing the QCD dynamics to determine the optimal 
smearing functions. 

We use interpolating fields which smear the 
static quark Q relative to the light quark q 



*.,(«,0 



J2siir)Qix + f,t)Tq{x,t) (1) 



and choose as a basis the set of N = (L/2 -+ 
l)(L/2 + 2)(L/2 + 3)/6 independent "S-wave" 
smearing functions 

^i(^) = 7^ y^ ^r.Rn (2) 
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7 , ^r,Rr, 

Reo 



(This generalizes for other representations of the 
cubic symmetry group O). For L^ = 20^ spatial 
sites, N = 286 is a manageably small number, 
and so the N x N correlation matrix 



a,(t) 



J2{o\^s,ix,mi^io,o)\o) 



(3) 



contains all of the information available from two- 
quark interpolating fields. Fixing the origin of the 
light quark propagator, taking advantage that the 
static-quark propagator depends only on a single 
spatial variable, and utilizing fast Fourier trans- 
forms allows the inexpensive computation of the 



